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The three primary ferroics, namely ferromagnets, ferroelectrics and ferroelastics exhibit corre- 
sponding large (or even giant) magnetocaloric, electrocaloric and elastocaloric effects when a phase 
transition is induced by the application of an appropriate external field. Recently the suite of pri- 
mary ferroics has been extended to include ferrotoroidic materials in which there is an ordering 
of toroidic moments in the form of magnetic vortex-like structures, examples being LiCo(P04)3 
and Ba2CoGe2C>7. In the present work we formulate the thermodynamics of ferrotoroidic materi- 
als. Within a Landau free energy framework we calculate the toroidocaloric effect by quantifying 
isothermal entropy change (or adiabatic temperature change) in the presence of an applied toroidic 
field when usual magnetization and polarization may also be present simultaneously. We also obtain 
a nonlinear Clausius-Clapeyron relation for phase coexistence. 

PACS numbers: 75.85.+t, 65.40.-b, 75.30.Sg 
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I. INTRODUCTION 

Ferroic phenomenon or the presence of switchable do- 
main walls in an external field is a direct consequence of 
a specific broken symmetry [l[ . Loss of spatial inversion 
symmetry results in ferroelectricity whereas loss of time 
reversal symmetry results in ferromagnetism. Ferroelas- 
ticity is a result of broken rotational symmetry although 
it remains invariant under both spatial inversion and time 
reversal symmetries. The fourth possibility corresponds 
to when both spatial inversion and time reversal sym- 
metries are simultaneously broken. This is the case for 
recently discovered ferrotoroidic materials [1J where the 
long-range order is related to an ordering of magnetic 
vortex-like structures characterized by a toroidal dipolar 
moment. It is important to mention that the ferrotoroidic 
order is also related to magnetoelectric behavior [3] which 
is one of the main attractions of multiferroics-materials 
that exhibit two or more ferroic orders simultaneously. 
This class includes magnetoelastics as well as magnetic 
shape memory alloys. 

In general, both electric and magnetic toroidal mo- 
ments can be defined within the context of electromag- 
netism 0, In the present paper we will exclusively 
consider the magnetic toroidal moment which is the only 
one that has the required spacial and time inversion sym- 
metries Q- We note that in recent years the study of 
the unique properties of ferroelectric materials carrying 
an electric toroidal moment has received considerable at- 
tention d 0| 

Ferrotoroidic domain walls have been observed in 
LiCo(P04)3 using second harmonic generation technique 
Q and a symmetry classification of such domains has 
been given recently [9j. Spontaneous toroidal moments 
have been attributed to exist in the multiferroic phase 
of Ba 2 CoGe 2 7 (BCG) E3 

and may be related to the 



observed unusual magnetoelectric effects [llj, |l2| . Single 
crystalline thin films of MnTi03 with an ilmenite struc- 
ture also exhibit a ferrotoroidic structure [l~3l |. Neutron 
polarimetry indicates that the magnetoelectric MnPS3 is 
a viable candidate for ferrotoroidicity |14| . The magnetic 
phase transition in BiFe03 implies the appearance of a 
toroidal moment (lo| Ab initio calculations suggest that 
the olivine material Li4MnFeCoNiP40i6 is possibly fer- 
rotoroidic [HI]. 

A sequence of possibly two ferrotoroidal phase transi- 
tions has been considered phenomenologically for Ni-Br 
boracites [I?) ■ Another physical realization of toroidal or- 
der has been considered as an interacting system of discs 
with a triangle of spins on each disc [18J. Both charge 
and spin currents can lead to a toroidal state [l9| . Recent 
observation of orbital currents in CuO through resonant 
x-ray diffraction provides a direct evidence of antifcrro- 
toroidic ordering |2Cf Ferrotoroidic materials exhibit lin- 
ear magnetoelectric effect; in fact, the toroidic moment 
is related to an antisymmetric component of the mag- 
netoelectric tensor {puj =fi Moreover, the toroidic 
moment can be viewed as a quantum geometric phase 
[2lj . Beyond the magnetoelectric applications toroidic 
materials can also act as novel metamaterials [22j . 

One way of realizing the physical consequences of 
toroidal moment is to characterize the thermal response 
of the material to an externally applied field that cou- 
ples to this moment and gives rise to measurable caloric 
effects. Actually, caloric effects are inherent to every ma- 
terial and are commonly quantified by the adiabatic tem- 
perature change or to the isothermal entropy change that 
occur when an external field is applied or removed. In 
the solid-state, the most studied caloric effect is the mag- 
netocaloric effect [2a |. mainly after the discovery in the 
mid-nineties of materials displaying giant magnetocaloric 
effect in the vicinity of room temperature 2J]. How- 
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ever, in recent years other caloric effects such as the elec- 
trocaloric |25j . elastocaloric (2(| or barocaloric [27[ have 
also received considerable attention. 

A crucial feature common to most materials exhibit- 
ing a giant caloric effect is the occurrence of a first-order 
phase transition. The expected large (discontinuous) 
change of the order parameter at the transition involves 
a large entropy content (associated with the latent heat) 
which is at the origin of the giant caloric effect. More- 
over, a strong coupling between different degrees of free- 
dom such as structural, magnetic, electric, etc. enables 
the transition to be driven by different fields conjugated 
to such properties. The study of caloric effects is thus a 
convenient method in order to study the thermodynam- 
ics of this class of complex materials. For instance, this 
should apply to multiferroic materials which are expected 
to show more than one caloric effect, e.g. simultaneous 
electrocaloric and magnetocaloric effects, or more pre- 
cisely a magnetoelectrocaloric effect. The latter has not 
been reported yet. The present paper deals with the ther- 
mal response and the associated caloric effects resulting 
from changes of toroidal order in ferrotoroidic materials 
(undergoing a paratoroidic to ferrotoroidic transition). 
These changes are likely to give rise to a toroidocaloric 
effect in this class of materials. Indeed, the study of 
toroidocaloric effect is expected to provide new insights 
into understanding the important problem of switching 
the toroidal moment 0- 

In this paper we study the thermodynamics of mul- 
tiferroics and first derive expressions for magnetoelec- 
trocaloric and toroidocaloric equations in Sec. II and Sec. 
Ill treating magnetization, polarization and toroidization 
as independent order parameters. In Sec. IV we present 
a Landau free energy for a paratoroidic to ferrotoroidic 
transition with a specific coupling between toroidization, 
polarization and magnetization, derive a phase coex- 
istence (Clausius-Clapeyron) relation and compute the 
toroidocaloric effect. Possible experimental implications 
are discussed in Sec. V in which we also present our main 
conclusions 



general Maxwell relations: 



/ dS 



f dXj 

\ dT 



(2) 



{Yj} 



Here {Xi} denote generalized displacements thermody- 
namically conjugated to the fields {Yi}. Any interplay 
between the different degrees of freedom will be taken 
into account through the fact that any Xi is, in princi- 
ple, a function of all fields. Therefore, the interplay must 
be introduced through the state equations. 

For each generalized displacement a caloric effect will 
occur when the corresponding conjugated field is varied. 
Suppose, for instance that the field Yi changes from 
to Yi so that the system passes from a state (T i5 0) to 
(T/,Yi), where T and T/ are the temperatures of the 
initial and final states respectively. The corresponding 
change of entropy is given by, S(Tf, Yi) — S(T i: 0). The 
two limits of interest that quantify the caloric effect asso- 
ciated to the property Xi conjugated to the field Yi, are 
the isothermal and the adiabatic limits. In the isothermal 
case Ti — Tf ~ T and the thermal response is character- 
ized by a change of entropy given by: 



A5(T,0->Y) = S(T,Yi) -S(T,0) 

^ /ax, 



V dT 



{Yj} 



dYi. (3) 



In the adiabatic limit, AS — 0, and the thermal response 
is quantified by the change of temperature given by: 



AT(T l; — s- Y) = T/CYO - H(0) 

/•Yi rp / 9X . 

~ ~J C \df 



dY^ (4) 



{Yj} 



Therefore, the caloric response of a material to a given 
field Yj will be given by {^r , ){-y y 



A. Examples 



II. CALORIC EFFECTS: GENERAL ASPECTS 

Consider a generic thermodynamic system and take 
temperature (T) and generalized forces or fields ({Yi}) 
as independent variables. Differential changes of the 
fields will yield a differential change of entropy [S = 
S(T, {Yi})] given by: 



c 

— ( 

T 



where C/T = (fy)| Y -} defines the heat capacity at con- 
stant fields {Yi} and we have taken into account the 



In the case of the magnetocaloric effect the correspond- 
ing expression for the isothermal entropy change is: 



A S (T,„^„, = f (|?) (h| .,H 

and for the adiabatic temperature change it is: 



(5) 



(6) 



where H is the magnetic field (strictly we should write 
/JoH = B) and M is the magnetization. Notice that 
by just replacing H by the electric field E and M by 
the polarization P the corresponding changes of entropy 
and temperature that quantify the electrocaloric effect 
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are obtained. Similarly, by replacing H by the stress a 
and M by the strain e we get the corresponding changes 
of entropy and temperature for the mechanocaloric ef- 
fect. Therefore, a barocaloric effect (a particular case of 
mechanocaloric effect) involving pressure p and change 
in volume AV^ is also expected. 



III. MATERIALS WITH TOROIDAL ORDER: 
THERMODYNAMICS 

The three basic moments of the electromagnetic held 
are electric, magnetic and toroidal moments. Therefore, 
we assume that the three ferroic properties are character- 
ized by the corresponding moments per unit volume, po- 
larization (P), magnetization (M) and toroidization (r) 
that can be assumed as independent (vector) order pa- 
rameters. Here the toroidization is assumed to originate 
only from the existence of magnetic toroidal moments. In 
fact, this appears to be the case for LiCo(P04) 3 Q and 
Ba2CoGe20 dHHl]. The corresponding thermodynam- 
ically conjugated fields will be the electric, E, magnetic, 
H, and toroidal, G, fields. Therefore, the "Thermody- 
namic Identity" for such a (closed) system reads: 



dU = TdS + E • dP + H • dM + G ■ dr, 



(7) 



where U is the internal energy per unit volume, T, tem- 
perature and S 1 , the entropy per unit volume. The field 
G coupling to the toroidization is related to the electric 
and magnetic fields through G — E x H. This choice, 
however, deserves some discussion. Note that the natural 
conjugated field of the magnetic toroidal moment is V x B 
Q. Nevertheless, since we are only considering homoge- 
neous macroscopic bodies in thermodynamic equilibrium, 
from symmetry considerations we assume that G is the 
appropriate macroscopic field that enables external con- 
trol of the toroidization. This is actually in agreement 
with references and [Ij where it is shown that the free 
energy of a system with magnetic toroidal moment must 
include a term proportional to the product r • G. It is 
worth pointing out that the field G cannot be modified 
independently of the fields conjugated to polarization and 
magnetization. From the Thermodynamic Identity (fT]). 
this relationship between G, E and H requires that, 



dU 
97 



S,P,M 



dU 
~dP 



S.M.T 



( — ) 

\dMJ 



(8) 



S,P,T 



We now define the Gibbs free energy Q through the fol- 
lowing Legendre transform: 



Q = U — TS EP-MH tG. 



(9) 



Differentiating this expression and replacing dU in ([5]), 
we obtain: 

dg = -5dT-[P + (Hxr)]-dE-[M+(rxE)].dH. (10) 



Taking into account that double differentiation of Q is 
independent of the order in which it is carried out, we 
obtain the following Maxwell relations: 

9e) th = (dr)~~ + UX (df, 



EH 



(11) 



and 



dS 



_ /dM\ 

T.E V ®T J h e 



dr 

dT 



EH 



x E. 



(12) 



H,E 



We now take into account that polarization and mag- 
netization can be written as the sum of an intrinsic term 
originating from (pre-existing) free electric and magnetic 
moments and a contribution arising from the toroidal 
moments. The toroidal contributions satisfy [U 13, Il9j : 

P t = -t x H, (13) 

and 

M,=txE. (14) 

Then, the Maxwell relations (jTTJ) and (fT2")l can be written 
in the form: 



dS_ 



T.H 



and 
dS 



dT 



fdMi 



EH 



dP t 

dT 



E,H 



dP 

df 



E,H 

(15) 



T.E 



V dT J h,e V dT 



fdM t \ _ 

H,E V dT / H ,E 

where the intrinsic contributions to the polarization and 
magnetization are related to the total as P = Pi + Pt 
and M = Mi + M t . 

A second set of Maxwell relations are obtained from 
differentiation with respect to E and M. This yields: 



fdPA _ f dM 



\dH ) 



T,E 



8E 



(17) 



T.H 



where we have taken into account that the intrinsic com- 
ponents Pi and Mi do not depend on H and E, respec- 
tively. This means that magnetoelectricity in the system 
originates only from the toroidal order. Notice that Eq. 
(jTTJ) just expresses that the magnetoelectric tensor a^- 
obeys Pi = \f ( /-.', • n, , //, and Mi = onjEj+XijHj, where 
Xij and xYj are dielectric and magnetic susceptibility ten- 
sors, respectively. From the r • G term in Eq. (9) it also 
follows that the components of the toroidal moment obey 
the proportionality: t\ ~ (a 2 3 - 032), i~ 2 ~ (a 3 i - ai 3 ) 
and r 3 — (ai2 — a%i). 

An interesting relationship between r and Pt, M t can 
be obtained considering that: 

P,xM( = -(r x H) x (r x E) 
= (t x E) x (t x H) 
= -(r x E) x (H x t) 

= —t X (E X H) X T 

= -txGxt. (18) 
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Taking into account the general vectorial relation A x 
BxA = A(A-B) — iyl 2 B, we can rewrite the preceding 
equation as: 



P* x M ( 



1 



t 2 G-t(G-t) 



(19) 



which points out that P t x M t is different from zero only 
under an applied toroidal field G. Here r = \t\ denotes 
the magnitude of toroidization. 



A. The toroido-caloric effect 

The equation quantifying the toroidocaloric effect un- 
der an applied electric or magnetic field is obtained using 
Maxwell relations and (TT21) . We obtain the following 
isothermal changes of entropy: 



A5(T, -> E) 



dT 



EH 



II 

dT 



EH 



dE (20) 



and 



AS(T, ^ H) 



H 



dT 



H,E 



Or 

df 



— — I x E 



H,E 



dH. (21) 



The corresponding equations for the adiabatic change of 
temperature are: 



AT(Ti, ->• 

rri 

C 



E) 



and 



AT(T,,0^ 
C 



dT 



H) 

'fdMt 



E,H 



E,H 



<m (22) 



V dT 



H,E 



9x\ 
9T/ 



— — I x E 



H,E 



• dH. (23) 



The second term in the square brackets in Eqs. (20) - (23) 
represents the toroidal contribution. It is interesting to 
notice that when applying an electric (magnetic) field, 
no toroidal contribution to the caloric effect will occur if 
the magnetic (electric) field is zero. 

If only toroidal order is present in the system the above 
equations for isothermal entropy change reduce to: 



AS(T,0^E) = — i 



d_Pt 

dT 



and 



AS(T,0->-H) 



H fdM t 



V dT 



E,H 



H,E 



dE 



dH. 



(24) 



(25) 



One finds similar expressions for the adiabatic tempera- 
ture change AT (T,-, — > E) and AT (T, H). 



IV. LANDAU MODEL FOR A MATERIAL 
UNDERGOING A FERROTOROIDIC 
TRANSITION 

We assume that the toroidization (r = |r|) is the or- 
der parameter to describe a paratoroidic to ferrotoroidic 
transition and propose the following minimal Landau 
model which includes a coupling between toroidization 
(t), magnetization (M) and polarization (P), and the 
presence of toroidal (G), magnetic (H) and electric (E) 
fields: 



F(T,t, P,M) 



iA (T-T c V + ^r 4 



2 Xp 



1 



~}M 2 



+ KT ■ (P x M) - G • T 

- H-M-E-P. (26) 

Note that the term k is the lowest order symmetry al- 
lowed term (satisfying space and time reversal symme- 
try) which provides the coupling between toroidization, 
polarization and magnetization (see [28^. Here the in- 
verse toroidic susceptibility Xt 1 = Aq(T — T®), where A 
is the "toroidic stiffness", T° is the transition tempera- 
ture and C > is the nonlinear toroidic coefficient. In 
particular, xt — dP/dH = <9M/<9E. We also note that 
a Landau theory of ferrotoroidic transitions in boracites 
[13 and Ba 2 CoGe 2 7 (BCG) [3 has been considered 
previously. However, these studies did not consider any 
caloric effects. 

Minimization with respect to polarization and magne- 
tization gives: 



and 



dF 
&P 

dF 
&M 



Xp l P - E + k(M xt) = (27) 



Xm 1 M-H + K (xxP)=0. (28) 



We now solve these two equations assuming (for simplic- 
ity) that E = (£,0,0) and H = (0,H,Q) and therefore 
G = (0, 0, EH) and r = (0, 0, r). For P [= (P, 0, 0)] and 
M [= (0, M,0)] we obtain: 



P = X p E-K Xp X m HT + 0(T 2 ) ~x P E-aH (29) 



and 



M = Xm H - KXpXmET + 0(r 2 ) ~ X mH - aE, (30) 

where in the above two equations we have neglected the 
nonlinear magnetoelectric effects. The magnetoelectric 
coefficient a = KXpXmT is a quadrilinear product of elec- 
tric susceptibility (x P = dP/dE), magnetic susceptibil- 
ity (Xm — <9M/<9H), the coupling constant k and the 
toroidization r. Thus, either for k = or r = there 
is no magnetoelectric effect. Substitution of P (f2"9"|) and 
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M (|30|) in the free energy ([2B1 gives the following general 
type of effective free energy: 



F e = F (E,H) 



1 * 
3 BT 



-Ct 4 + At, 
4 



with 



F 
A 

T c 

D 
A 



1 



A Q (T -T°) - K 2 XpXm [x m H 2 
A (T-T C ), 

T° + ^-XpXmlXmH 2 + X pE 2 } 

3K 3 x 2 mX 2 p EH, 
(KXmXp ~ l)EH. 



XpE 2 



(31) 

(32) 

(33) 

(34) 

(35) 
(36) 



This corresponds to the free energy of a system sub- 
jected to an effective external field A, proportional to 
the toroidal field G = EH. When G = 0, and there- 
fore from (|35"1) B = 0, the free energy (|3Tj) describes 
a paratoroidal-to-ferrotoroidal second-order phase tran- 
sition whereas under the application of a toroidal field 
G ^ (and B^O), the transition becomes a first-order 
one. Actually, the physics contained in (l31"j) is very rich 
since, in addition to A, the cubic B coefficient also de- 
pends on the toroidal field G and the linear A(T) term 
explicitly depends on the E and H fields. In other words, 
A, A(T) and B are not independent. This leads to a com- 
petition between B and A depending on the value of the 
coupling constant k (i.e. the choice of the ferrotoroidic 
material) and to a nonlinear Clausius-Clapeyron equa- 
tion. It is worth pointing out that the addition of non- 
linear terms in eqs. |29|) and (|30|) would lead to higher 
order terms in the expansion (I3ip that go beyond the 
minimal model. However, within the spirit of the Lan- 
dau Theory, we expect such terms are not essential. 

For some purposes, it may be convenient to rescale the 
free energy according to the definitions proposed in [2!| : 



«o 
h 

fr 



A C 
B 2 ' 
AC 2 



|S|3 

c 

W\ T 

F eff 

B 4 



C 3 



(37) 



which yields the following rescaled free energy, valid for 
B^O only: 



1 2 1 , 1 

fr = fo + g a T r + 3 T r + J T r 



hr r 



(38) 



As usual, the liner coefficient a is temperature dependent 
a = ao(T — T c ) and provides the temperature scale, while 
h provides the scale of the external effective field. 



A. Phase diagram 

Figure Q] shows the phase diagram for the rescaled vari- 
ables (|37p defined above. In the upper panel we have 
plotted the behavior of h vs. a(T) which gives the co- 
existence line. The region above the line is paratoroidic 
whereas below it is ferrotoroidic. The coexistence line 
ends in a critical point that can be calculated from the 
condition: 



dh 



= a + 2t + 3t 2 = 0. 



(39) 



One obtains 

(/i c ,Q c ) (27:3) 

k* < 1 whereas for k 



the critical point is located at 
Notice that the field is h < for 
> 1 we have h > but 



h < h c = A . The existence of this critical point is also 
revealed from the turning point in the behavior of r r 
vs. h shown in the lower panel. Below the critical field 
(h < h c ), there are two possible values of t t related 
to the two possible wells in the free energy, as it is 
schematically illustrated in the insets. Different symbols 
denote that the results have been obtained for different 
values of the effective coupling parameter n* = n X mXp 
by solving numerically the model (|3 1[) . The subsequent 
application of scaling relations defined in (13"T1) leads to 
the curves shown in Fig. [TJ 
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FIG. 1: (color online). Rescaled phase diagram as a func- 
tion of the field h for the temperature coefficient a(T) (upper 
panel) and the order parameter tv (lower panel). Different 
symbols denote results obtained for different values of the ef- 
fective coupling constant k* = nXpXm- 
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Actually, the unsealed free-energy (13TT) is more suit- 
able for studying independently the effect of the external 
toroidal field G and the material dependent coupling pa- 
rameter re. For convenience, we set Xp = 1j Xm = 1; 
C = 1, = 1 and Aq = 1, which renders the following 
dependence for the model parameters: 

S = 3k 3 EH = 3re 3 G, (40) 

A = (k - l)EH = (k - 1)G, (41) 

A(T) = T-T c , (42) 

T c = 1 + k 2 (E 2 + H 2 ). (43) 

The only free parameter is the coupling constant re. From 
now on we will restrict to values of k > and conse- 
quently to B > 0. Under these conditions, different situ- 
ations can be considered: 

• B > and A = 0, corresponds to n = 1. In the low 
temperature regime the solution corresponds to a 
minimum located at r < 0. 

• B > and A > 0, corresponds to n > 1. In this 
case the low temperature minimum also occurs at 
t < 0. 

• B > and A < 0, corresponds to < n < 1. In that 
case, competition between r < and t > occurs: 
A < favors the minimum to occur at r > while 
B > favors the minimum to occur at r < 0. 

As a reference case, we also consider the possibility of 
G = 0. In that situation, 

• £ = and A = for all values of n. The tran- 
sition is continuous and occurs at T£ = 1. In the 
low temperature phase the free energy shows two 
symmetric minima at ±tq. 

The corresponding phase diagram obtained from the un- 
sealed free energy ([31]) is depicted in Fig. [2j It shows 
the behavior of the effective field A as a function of the 
transition temperature T t (or A(T t ) = T t — T c ) for the 
three different regimes of the coupling parameter men- 
tioned above; namely re > 1, re = 1 and k < 1. (i) For 
k > 1 (above) the transition exists for values of A above 
the critical field A c (re) > 0, which in turn depends on 
the value of re. Moreover, T t increases with A as it can be 
seen from the representative curves obtained for re = 1.05 
and 1.10. (ii) For k — 1 the effective field is A = and 
the coexistence line is horizontal starting at the critical 
point (at the origin) corresponding to the limiting case 
of G — (B — 0) and denoted by a blue dot. (iii) Fi- 
nally for k < 1, one has A < and the transition exists 
for values of A arbitrarily close to zero at temperatures 
T t < T c . For increasing values of G, A decreases while 
A(T t ) increases; thus eventually the latter may become 
positive. We stress that all curves correspond to discon- 
tinuous phase transitions except the point located at the 
origin that corresponds to B — and therefore to a con- 
tinuous phase transition. 




A(T t ) 



FIG. 2: (color online). Phase diagram showing the effective 
field A as a function of the transition temperature A(Tt) = 
Tt — T c for representative values of « > 1, K = 1 and re < 1. 
The inset shows the coexistence line in rescaled variables. The 
(blue) dot on the A = line corresponds to the case of B — 0. 



Figure [3] shows the temperature behavior of the 
toroidal order parameter r for three representative val- 
ues of the coupling parameter namely k = 0.90, 1.0, 1.05 
and different values of the toroidal field G = EH as in- 
dicated in the lateral panels. The case k = 1.05 (lower 
panel) nicely illustrates the effect of the field on the fer- 
rotoroidal transition in the region of K > 1. For zero-field 
(G = 0), the transition is continuous and, indeed, occurs 
at T = T° = 1 (symmetric curve with double branch). 
By increasing the field (G > 0), the transition first disap- 
pears (continuous cross-over, no singularity, from r 
to r = 0) and subsequently for higher values of the field 
a first-order transition occurs at a temperature T t > T® 
revealed by a discontinuous jump in r (indicated by an 
arrow), whose magnitude increases with G. In the case of 
K = 1 (middle panel), the transition is discontinuous and 
exists for every value of G ^ 0, although it is continuous 
for G = 0. In the case re < 1 as soon as the field is ap- 
plied a first-order transition from r < to r > occurs 
with cooling. For low values of G, the transition takes 
place at T t < T c but with increasing values of G, T t also 
increases and the transition takes places above T c . This 
behavior results from the competition between B and A. 

The temperature behavior of the polarization P and 
the magnetization M can then be obtained from the ex- 
pressions (f2"9"]) and (|30p . respectively. The results are 
shown in Fig. 0] for the case of re = 1.0 and different 
values of G = EH (E = 1.0). 
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FIG. 3: (color online). Toroidal order parameter as a function 
of temperature for k — 0.90, « = 1 and k = 1.05 and selected 
values of the applied toroidal field G. The arrow in the lower 
panel indicates the occurrence of the discontinuous transition. 



FIG. 4: (color online). Toroidal order parameter (r), magne- 
tization (M) and polarization (P) as a function of tempera- 
ture A(T) = T-T c for n = 1.0. The toroidal field is modified 
by changing the magnetic field H but keeping constant the 
value of the electric field E = 1.0. 



B. Nonlinear Clausius-Clapeyron Equation 

Since the transition predicted by the model is discon- 
tinuous, it can be characterized by means of the corre- 
sponding Clausius-Clapeyron equation. Such equation 
relates the slope of the coexistence curve (Fig. [2j to 
the magnitude of the discontinuities in the order param- 
eter and entropy at the transition temperature. We re- 
call that in the present model, the harmonic (A(T)) and 
cubic (B) coefficients and the effective field (A) are not 
independent. As noted below, this will give rise to a 
Clausius-Clapeyron equation which turns out to be non- 
linear in the order parameter discontinuity. 

Indeed, at the transition, one can write: 



partial derivatives are: 



dFA „ A Q _ 2 



dF e 



ox 



E(k-I) \dH J T ' 



(46) 
(47) 



where S is the entropy and the derivative with respect to 
the field A is taken at constant E. 



V OH 



= -H(l + A k 2 t 2 ) + E (k 3 t 3 + t(k - 1)) 



(48) 

Taking into account ([35]). (jT7j) and (gSJl, condition (|33|) 
reads: 



dF, ' « = dF ' 



(44) 



where (I) and (II) are the para- and ferrotoroidal phases 
respectively. For each phase, the derivative of the effec- 
tive free energy can be written as: 



dF„. = 



dF e 
dT 



dT 



dFe 



dX, 



(45) 



where F e (T, A) is the thermodynamic free energy given 
by the expression (|31[) with r = t(T) in equilibrium. The 



dT(Su-Si) = d\(AF e )dX 
dX 



with 



E(k - 1) 



\anJ T 



\dH] T \dHj T 
+ (r(/Z) -T(j))(K- l)E. 



(ii) 



\ dH J T 



ii)) H 



(49) 



(i) 



(50) 
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FIG. 5: (color online). Clausius-Clapeyron equation. We 
have plotted the coexistence line (black dots) for n — 1.10 
obtained at constant E = 1.0. Squares in red denote the 
results of computing the right hand side of Eq. (50) from the 
discontinuities of r and S at the transition point. Finally, 
the continuous line (blue) is the numrical derivative of the 
coexistence line. 



Finally, the corresponding Clausius-Clapeyron equa- 
tion takes the form: 



applied toroidal field (G — EH) has been plotted in Fig. 
[5] for the same three representative values of k as in Fig. 
El As it can be observed, with increasing field value the 
jump in AS increases. 




1.0 1.5 
T (reduced units) 



dT 



AS 



dx(AF e 



(S H - S!)E(k - 1) 
V dH It 



\dH IT 



(51) 



Notice that for k = 1, the slope of the coexistence curve 
is zero although there exists a jump in both r and S at 
the transition point. In Fig. [5] we have plotted the coex- 
istence line for k = 1.10 (black dots) and its numerical 
derivative (continuous blue line). The squares (full red) 
are the results of computing the right hand side of Eq. 
(|5ip from the magnitude of the discontinuity of t and 
S at the transition temperature. The agreement is very 
good. 



C. Toroidocaloric effect 

The toroidocaloric effect can be computed as: 

f G=EH / Q • 

S(T, G = EH) - S(T, G = 0) = / I — I dG 



dT 



G 



1 



= ^A [t 2 (T,G = 0)-t 2 (T,G = EH)]. (52) 

This is an important equation for the isothermal toroidal 
field-induced entropy change which is analogous to the 
expression found long ago for the adiabatic electric field- 
induced temperature change in the case of the elec- 
trocaloric effect [3(3]. The isothermal entropy change 
(AS) as a function of temperature for various values of 



FIG. 6: (color online). Toroidocaloric effect as a function 
of temperature and selected values of the maximum applied 
toroidal field, k = 0.90, 1.0 and 1.05. 



It is interesting to relate this entropy change charac- 
terizing the toroidocaloric effect with the corresponding 
changes giving the magnetocaloric and electrocaloric ef- 
fects. This can be done using Eqs. (|2"9")) and (|3"0")) , a 
straightforward calculation gives: 

AS(T, -> G = EH) 

' - [AS P (T, H,0^E) + AS M (T,E,0^H)} 



^XpXm 



(53) 



where ASp and ASm are the entropy changes giving the 
electrocaloric and the magnetocaloric effects at constant 
magnetic and electric fields, respectively. 

To apply our results to a specific material we need to 
obtain parameters (e.g. k) for LiCo(P04)3 using the data 
from second harmonic generation [8j and other thermo- 
dynamic measurements, in particular the toroidic sus- 
ceptibility. The latter is related to the antisymmetric 
part of the magnetoelectric tensor as discussed in [lj| for 
BiFeOa. It is interesting to point out that Eq. ([53]) sug- 
gests that if the electric and magnetic susceptibilities are 
known, k could be estimated from measurements of the 
entropy changes AS, ASp and ASm- 
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V. CONCLUSIONS 

With the discovery of the fourth kind of primary ferroic 
materials, namely the ferrotoroidics 0,3, it is important 
to understand the equilibrium thermodynamic proper- 
ties of such materials. We provided a basic framework 
to calculate the toroidocaloric effect based on a Landau 
free energy. Our main finding is the isothermal change 
in entropy as a function of the applied toroidal field for 
different values of coupling between toroidization, polar- 
ization and magnetization. The fact that the applica- 
tion of a toroidal field modifies the order of the tran- 
sition from continuous to first-order is very informative 
regarding caloric effects since larger changes of entropy 
are expected to be induced thanks to the large entropy 
content associated with the latent heat of a first-order 
transition. However, when dealing with real materials, 
the existence of energy losses arising from hysteresis and 
domain wall effects [23|, Hl[ (which are intrinsic to first- 
order transition and are expected to reduce the magne- 
tocaloric efficiency) should be taken into account. This 
important aspect has not been considered in the present 
paper since strict equilibrium situations are assumed. In 
any case, our predictions should be observable in exper- 
iments on materials such as LiCo(Po4)3, Ba2CoGe20y 
(BCG), MnTi0 3 , MnPS 3 and some boracites. Below 21.8 
K neutron diffraction measurements in LiCo(Po4)3 indi- 
cate simultaneous presence of ferrotoroidic and antiferro- 
magnetic (AFM) order as a result of Co 2+ ion ordering 



hibits ferrotoroidic behavior above 6 T magnetic field due 
to spin flopping fl3l ] . These results confirm the existence 
of a ferrotoroidic phase in a number of materials, how- 
ever, the measurement of the toroidal moment as a func- 
tion of temperature and toroidal field has not been under- 
taken. Therefore, at the present stage it is not possible to 
contrast our predictions with experimental results. We 
expect, however, that our results including the nonlin- 
ear Clausius-Clapeyron relation will be a motivation for 
experimentalist to undertake experiments aimed at char- 
acterizing the thermodynamic behavior of ferrotoroidic 
materials in the vicinity of the phase transition. 

Using the Landau free energy we can also obtain 
the profiles of ferrotoroidic domain walls by including 
symmetry allowed gradient term (Vt) 2 [19j . i.e. the 
Ginzburg term. Such domain walls have been observed 
in LiCo(Po4_)3 using optical second harmonic generation 
techniques [8] . With doping induced disorder in such ma- 
terials we expect that novel phases such as toroidic tweed 
and toroidic glass should also exist and remain to be ob- 
served experimentally. With symmetry allowed coupling 
of strain to toroidization, if we apply stress to such a 
crystal we expect toroidoelastic effects, i.e. a change in 
toroidization with hydrostatic pressure or shear. These 
important topics will be explored in near future. 
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